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ABSTRACT 


Combat models often involve target detection times which may vary with 
different observers due to characteristics of personnel, or detection systems. They 
may also be affected by different environmental factors such as visual levels, sea 
states, terrains, etc. There is often interest in quantifying the effects of different 
observer characteristics and environmental factors on detection times. A hierarchical 
gamma/Weibull regression model is considered which can incorporate observer 
characteristics and environmental effects which may influence the time to detect 
targets. Numerical procedures for the estimation of parameters of the hierarchical 
‘ gamma/Weibull model based on maximum likelihood are described. Results of 


simulation experiments to study small sample behavior of the estimates are reported. 
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I. INTRODUCTION 


A. BACKGROUND: 


Combat .models often involve target detection times which may vary with 
” different. observers due to characteristics. of personnel, or detection systems. They 
may also be affected by different environmental factors such as Visual levels, sea 
States, terrains, etc. There is often interest in quantifying the effects of different 
observer characteristics and environmental factors on detection times. In Chapter 
‘TI a hierarchical gamma/Weibull regression model [Ref. 1] is considered which can 
incorporate observer characteristics and environmental effects which influence the 
time to detect targets. Numerical procedures for the estimation of parameters of the 
hierarchical gamma/Weibull model based on maximum likelihood are described. 
Results of simulation experiments to study small sample behavior of the estimators 
are reported. In the remainder of this chapter, two models related to the one 
considered in this thesis will be described. Numerical procedures for parametric 
estimation for these models are studied in Bae [Ref. 2]. Experience gained 


concerning the numerical stability of these procedures will be described. 


B. TWO RELATED MODELS 
Two parametric models for the distribution of time to detect targets were 
considered in a thesis written by a Korean officer,Mr.Bae in Sept.,1989 [Ref. 2]. 
Mode! 1 :[Ref.2:pp.2-14] supposes that there are M observers. Observer i has 
parameter A, which characterizes the ability of the observer to detect a target. 
Observer i is presented with N; targets. The parameters of all M_ observers, 


Aj,A>,...Ay are assumed independent with a common two parameter gamma 


‘distribution. Given A, the detection times of observer i afe conditionally 
. independent, having Weibull distributions with known shape parameter e= and 
known scale parameter Mi for j=1,2,.,N; Some of the detection times may be 
censored. Bae [Ref. 2:pp. 2-20] reports on numerical procedures to find the 
maximum likelihood estimates for the shape and scale parameters of the gamma 
distribution. 
Model.2 : is a Weibull regression model [Ref.2:pp.21-30]. It assumes observer 
i has explanatory variables x;;, Xj2, .... Xjp relating to his j" target representing 
factors which influence his time to detection. The detection times for the observers 
are independent random variables having Weibull distributions, Again, some of the 
detection times may be censored. The scale parameter of the Weibull distribution for 
the detection times of the j" target by observer i is of the form per Vv where 
xf-D 390, ; the shape parameter is of the form e . Bae [Ref. 2:pp. 21-36] 
reports on numerical procedures to find the maximum likelihood estimates of {é,} 


and {B,}. 


C. OBJECTIVE ANP METHODOLOGY IN THE COMBINED MODEL 

The main effort of this thesis will be to study a full hierarchical gamma/Weibull 
regression model which is a combination of the two previous models, Since the 
models considered in Bae [Ref. 2] will be special cases of the one considered in this 
thesis, it is expected that instabilities found in the numerical procedures to estimate 
the parameters in the two previous models will appear for the combined model. Two 
numerical instabilities that result from the numerical procedures of Bae [Ref. 2] are 
described below. 

In Model 1, a modified Newton-Raphson method [Ref. 3] is applied to solve 


the nonlinear system of equations for the maximum likelihood estimates of the 


gamma ‘parameters. It is found that, particularly for small sample sizes, the slope of 
the likelihood surface is.very flat resulting in. numerical instability for estimation of 
the eae ‘parameter-y. As a result, in the full hierarchical model presented here we 
have used a single parameter gamma distribution as the second stage model. The 
maximum likelihood estimate of the parameter is found by the bisection root search 
method: [Ref. 4]. - 

A numerical overflow phenomenon was occasionally found in the Weibull 
regression model while iterating the Newton procedure for {€} in Model 2 
[Ref.2:p.25]. Once again, the flatness of the relationship determining the estimate 
of {€;} resulted-in this instability. As a result, the Newton procedure to determine 
{é} in the gamma/Weibull regression model of this thesis has some checks to detect 
this numerically instability, 

In Chapter II, a numerical procedure is presented to find the maximum 
likelihood estimate for a single parameter gamma version of Model 1. In Chapter III, 
the full hierarchical model is presented and numerical procedures to estimate the 
parameters given. In both chapters, simulation results are presented to study the 
small sample behavior of the estimates. All simulations were carried out on an IBM 
3179 G mainframe computer at the Naval Postgraduate School using the APL 
GRAFSTAT random number package [Ref. 5],[Ref. 6]. 

We hope that the hierarchical gamma/Weibull regression model and these 
estimation procedures will be a useful tool to describe and predict target detection 


times, which are one aspect of the effect of human performance on the battlefield. 





11.4 SINGLE PARAMETER GAMMA HIERARCHICAL MODEL 


A, MODEL. DESCRIPTION 

This model ‘is very similar to Model’ 1 in [Ref. 2]'. The difference is that the 
parameters A,, i = 1,2,....M, which reflect the abilities of observer i to detect targets, 
are assumed to be independent having a single-parameter gamma distribution, 
GAM(o,a) rather than having two parameters. The form of gamma density function 
used here is 


(a8) 


-a8 2.1 
Ta) e for §6>0. (2.1) 


g(8) = 


For numerical reasons we will parametrize a as e". For convenience we will 
reiterate the remainder of the assumptions of Model 1 given in [Ref. 2]. Given A, 
= 6, the times it takes for observer i to detect target j, J = 1,2,...,N;, denoted by Uj 
are assumed to be conditionally independent random variables with Weibull 


distributions , WEN p,,e°* having cumulative distribution function 


P(U,st |A,=6}=1 - exp{-e(d/u,)* } t20 (2.2) 


independent of other observers. The variations of the A, are introduced to represent 


the individual differences between the observers. When the i" observer is presented 


‘with his j" target, only an opportunity time O, is allowed for him to detect it. An 


observer either successfully detects the target within this time or never detects it. 


1 In the model of [Ref. 2], a two parameter gamma distribution is used. 








Data. for the i observer consist of times of detection for the successes and the 


lengths of opportunity times-for-the failures. For each i = 1,2,..,M,j = 1,2,...,N; , 


let 
Y, = min(fnU,,InO,\ (2.3) 

and 
_J 1 #U,sO, 2.4 
ey -{ 0 otherwise ° a 


The Y;; are the censored In-detection times and A, is an indicator of whether the In- 


time to detect the j" target by the i observer is censored or not. Let 
J get by 


c.= DA, (2.5) 


be the number of targets detected by observer i. In this chapter, we will assume the 
Weibull parameters {u,} and {¢,} are known constants, and we are only interested 
in estimating the gamma parameter 7. In the next section, the likelihood function 
for n is given and the bisection method for finding the maximum likelihood estimate 


of n is described. In the final section, we present results of a simulation study of the 


behavior of the estimators for small sample sizes. 





B. MAXIMUM LIKELIHOOD ESTIMATION AND THE BISECTION ROOT 
SEARCH 
4. The Likelihood Function 
Given.A, = 8, the conditional likelihood function [Ref.2:p.3] for observer 


i using the censored In-times Yj is 


P _ 
L(wEle) = T[lecOr**e}exg-ectrtr"], 6) 
jal 
Let 
N; % 
S, = 2 exP\0y - inp, je} (2.7) 
j= 
and 
N 
K, = en EAilyy-nn,)e™ - é|| , (2.8) 
Phd 


Equation (2.6) can be rewritten as follows : 


L{u.E |0) = 6° Kexp(-05) . (2.9) 


From the equations (2.1) and (2.9), it follows that the unconditional likelihood 


function for observer i is 





Afurkad) = [{L{u.&|\e(@ae 


Be oC) ae 
fre%e aoe (2.10). 
2 ork 
= Sey (a+k) . 


Ye 
(Si, 00)" b= 


‘C1 
In equation (2.10), if C, = 0; then JJ («+k) = 1 . Recall the parametrization o 
20 


=e", The unconditional In-likelihood function for observer i can be rewritten as 


InL(y.t.n) = 
C;-1 (2,11) 
InK, + e"yn - (C,+e")In(S,+e") + LU In(e"+hk) 
k=0 


C1 
where if C; = 0, then > In(e"+k) = 0 . Since the observers are independent, the 


k=0 
unconditional In-likelihood for all M observers is 


re 
i (2.12) 


i=l 
M C1 
> {nur n-(C,+e") In(S,+e")+ 3 nero ‘ 
k=0 


In this chapter, we will focus attention only on finding the maximum likelihood 


estimate of n. The derivative of the In-likelihood with respect to n is 














= fq) * 
on ? 
ut 5-¢: . C1 ‘ (2:13) 
e"y {mn + = — - ln(S;+e") + 
i=l | S,+e" kd e" +k 


The maximum likelihood estimate for n is the-solution of the equation f(n) = 0; 
| 2. Bisection Root Search for 7 

Initially, the Newton procedure was used to numerically solve the equation 
f(n) = 0, Unfortunately, the procedure frequently either convergéd to an 
unreasonable number or encountered problems of numerical instability. Plots of f(n) 
indicate that f can have multiple zeros, A:representative plot of f is shown in Figure 
A ? (next page). In Figure A, one v symbol indicates a zero of f, while two ¥v 


symbols indicate the reasonable root to be used as the estimate for n. 


2 This is a typical graph for f{(n) which was generated for simulated data with O=10, M=15, N=15 
and random seed = 16807. The true value of ») is 1. 
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Figure A Multiple Roots of the Derivative of the In-Likelihood Function for 


ch method [Ref.4] 
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As a result of these possible mult 


| the estimation 


ai 


The following sub-section will det 


=0. 


d to solve f(n) 


1S use 


procedure for n. 





>: a nitialEstimate forn 
The initial.estimate for n is-similar to its initial estimate-in the two- 
parametér gamma model [Ref. 2], and is obtained by setting 3B, = 0 in that 


" procedure. Consequentially,.the initial estimaté of n is 


n° = “fiat ' (2.14) 


where 


th, = -0.5772 


MN, sea 
LEA loy-me)e | (2.15) 
aA — s j= 
Mm, rs MN eS 
DLA, 
iw ljut 


If ar <0 ,setn? = 04, 
b. Determination of Rounds for the Root 
The purpose. of this step is to bound the root that will be used as the 
estimate between two points, called n, and n,,. One of these points is n° and the 


other is found as follows : 


3 In MrBae's model [Ref2J, the regular gamma distribution is parametrized with shape parameter v=e"**e 


and scale parameter a=e". Setting 6)=0, a=v=e" yields the one-parameter gamma distribution. 
P g y Pp 


2 
. ~ aA v9 . oye . * 
4 The estimate titi“ <0 suggests that the variability of the data is small, Hence, it 


might be better to set 4° equa! to a large positive number, such as 13. A topic for future research is 
to explore better initial values for in this default case. 
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Hy=nookxd 


Set i sfilekd.  ke103..... 


where: is:a-‘small increment. 
Compare the signs of f(7,): and £(np) with f(n°), until the first different sign from f(n°) 
is.found. This yields a bound.on the closest root. to-n°. 
c. Bisection Search 
After determining the left bound, n, and right bound ,np, we begin 

the bisection search by picking the midpoint between them, called'7,,, and computing 
the value of f(ny). If £(m,) and f(n,). have the same signs, the root is on the other 
side of ny; thus, move the left bound n, to midpoint ny. If the signs are different, 
move the right bound to the midpoint, Iterate the bisection search until 
7Tr 

"Mm 

C. IMPLEMENTATION OF THE SIMULATION AND RESULTS 


<0.01 . The maximum likelihood estimate of 7 is ny. 








1, Simulation 
In order to study the sampling properties of the maximum likelihood 
estimate of n, different numbers of observers (M), targets (N,), and opportunity times 
(O;) are used as inputs to the simulation. Each set of inputs (M,N,,O,) has a 
simulation run .of (R) 100 replications and each run begins with the same initial 
random seed (466801743). Statistics, of mean bias (M.B), mean square error (M.S.E) 
and their standard errors ( S.E(M.B) and S.E(M.S.E) ),.are computed according to 


following definitions : 


R 
MB = 2S (in) (2.16a) 


ral 


it 



















Sate te oo : — | 3 
SELB) = |= (4,-1)-MBP (2:16b) 
NY RR-2) 


% R | 
pe oa -_ M.S.E = iG ony 278) 


rel 





fee cn area (2.17b) 
Rep n)-M.S.E) 


jc ; S.E(MIS.E) = . 


_ where 4, represents the point estimate of parameter 7 in the r" replication. The 
fraction: of detection level or the averaged uncensoring level (UC) for R (100) 


replications is determined by 


an (2. ) 


where C;,(r) is the number of targets detected by observer i in the r" replication and 
N,(r) is the number of targets presented to observer i in the r replication. 
Theoretically, the longer the opportunity times, the higher the UC level should be. 


An outline of a replication in the simulation is as follows : 
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. Input tO. Simulation 


M the number of observers (Possible values are 15 , 25 and 35 ) 


N, the number of targets. presented tothe ith observer , where i = 
1,2....M_ (Possible values are 15 ,25 and.35) 


; the SPpOREHTY time for ith observer to detect the jth 
target,(Possible values are 10 and 15). All of the opportunity 
times are the same. 


§; _ the shape parameter of the Weibull distribution for ith observer, 
where 1 = 1,2,..M ; ( Set €; = 0 for alli) 


ui; the scale parameter of the Weibull distribution ,WEI(u, , e*); 
(Set eal = 42 ) 


n the parameter for the single-parameter gamma distribution , 
‘GAM(a,a) , with a = e"; (Set the true value of n = 1) 


ion ri server 


Generate independent single-parameter gamma rv.’s A; from 
GAM(a,a),having the density function as in equation (2.1). 


Generate independent exponential distribution r.v.’s W, with mean 
1 ive, Wi; ~ EXP(1) for all i and j 

W, ef 
Compute the recorded In-times Y;; as in equation (2.3). 


Compute the target detection times: U, : 


Compute A, as in equation (2.4). 


W,, efi t » ee 
— | st|A=6) = P W,,<8) — = 1-exp{-6] — which is the 
4; Bi By 


conditional Weibull distribution of equation (2.2). 
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c. Initial n-and‘Bisection Search 
¢ Compute mm, asin equation (2,15).and determine the initial value 
of nas in.equation (2.14). 


¢ Perform the Bisection root ‘search: procedure for the maximum 
likelihood estimate. of 7, shown as the flowchart in Figure B. 


All these procedures were written.in APL codes,named "SIMULA1", 
"GAMMA", and "FVALUE"’, which are listed in Appendix C. 














“KEK 41 
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Figure B The Algorithm of Bisection Search 














2. Results. 

The results of the simulation experiments appear in Tables and graphs in 
Appendix A. Tables 1 and 2 contain the general statistics for the bias of | at the 
opportunity times O = 10.and 15 respectively. Table 3 in combination with box 
plots, shown in Figures 1 through 3, presents statistics of the bias of the estimates of 
nas a function of the numbers of targets and observers in the case O = 10, Figures 
4 to 6 show histograms of the bias with the number of. targets fixed and the number 
of observers changing. Similarly, Table 4 and Figures 7 through 12 focus on the case 


O = 15, In summary, the simulations indicate the following : 


° oy a time (O) 
ables 1 and 2 indicate‘that a longer Spporatty time is-reflected in 
a higher uncensoring level (UC), which slightly decreases the mean 
bias (M.B) and mean square error (M.S.E) of n. 


e The Number of Observers (M) 

Increasing the number of observers for a fixed number of targets 
results in greater changes in the mean bias and mean square error of 
n than increasing the number of targets with a fixed number of 
observers, All the box plots and histograms display this tendency. 
This phenomenon is quite reasonable, because the parameter n in 
gamma distribution reflects the variability in abilities of observers. 
Hence data with more observers will provide better estimation of n 
( less bias and standard error ). 


e The Number of Targets (N;) 

The box plots indicate that increasing the number of targets for a 
fixed number of observers has a small effect on M.B and M.S.E of n 
with no systematic trend. 











JI, THE HIERARCHICAL -GAMMA/WEIBULL REGRESSION. MODEL 


A. MODEL DESCRIPTION 

Recalling the assumptions of the model, we suppose there are M observers, 
indexed by i, The i” observer is presented with N, targets, indexed by j. The 
variable Uj; is the time it takes for observer i to detect the target j. Let Xj: , Xjo seo 
Xj) be the values of environmental variables (¢.g. terrain, sea state , atmospheric 
condition etc.) which may affect Uj, Given the quantified value of the ability for 
observer i, A; = 6, the U;’s are assumed conditionally independent random variables 
having Weibull distributions, WEI, ,¢°*+) , where the scale parameter, instead-of 


‘being constant, has the form 


By = exp{x,B} = eof rah. on) 


The regression coefficient 8, is the contribution to the Weibull scale parameter of 
the k'" explanatory variable. £, is usually the constant term in the regression. The 
random variables {A,} are assumed independent from a single parameter gamma 


distribution, GAM(a,a) with a = e", 


B. MAXIMUM LIKELIHOOD ESTIMATES FOR ALL PARAMETERS 
The model assumptions detailed above are the same as those in Chapter II, 


except that the scale parameter p, in the Weibull distribution is now a function of 
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- 
a 
“ 


the explanatory. variables. In this section we consider the estimation of all the 
parameters of n , {6;} and’{é}.° : 
From equations (2.7), (28), (2.12) and equation (3.1), we can rewrite the. 


unconditional In-likelihood function for all M observers .as follows : 


InL = 

M 5 C-1 
g(n.£,B). 5 in + en - (C+e")In(S.+e%) + Lin(e+h | 
BIS ke0 


ek | 


of tM ‘ (3:2) 
Where K, ~ ex} EAyy-x, Be i. é, | 


N, 
“ = “Gr \. 
and, aoe {oy x Be \; 


Yj, 4, and C; are defined in equations (2,3), (2.4) and (2.5) respectively. 

In the following sub-sections, we will describe the procedures which are 
employed to estimate all the model parameters. The general procedure includes the 
initial rough estimation of {8,} and {€} without the hierarchical gamma rwv.’s; a 
variability check for the need to include the hierarchical gamma parameter n; and 
finally, a large recursive procedure for estimating all the model parameters if a 
hierarchical model is necessary. 

1, Initial Estimation of a Simple Weibull Regression Model 

An initial model for the data is the simple Weibull regression model for 


the target detection times described previously as the Model 2 of Bae [Ref. 2:pp. 21- 





30]. The key steps of the estimation procedure are summarized as follows : 


§ In the Model 2 of Bae (Ref. 2], no gamma variable was involved. Only the Weibull parameters 
{B,} and {€;} were considered. 


18 


Step-LL : Initially, set? x,8°= InUy and €/ = 0. 
Step-1.2 : This is an.iteratively re-weighted regression step; cf. McCullagh 





and .Neélder [Ref. 7]. 
. Let 
r= ba yA" ee 
Wy = expiry} : (3.4) 
yy = Wye where h=0,1,..,.p (3.5) 
Zz A; wi] 5 ynPa & (-A,+ i) ew ety ; (3.6) 
Wy Wy v 


Regress the dependent variable z, against the independent variables uj, 


for h = 0,1,...,p ; e.g. the estimate has the form 


B= (UUWZ (3.7) 


where 


7 From this page on, all parameters with superscript 0 imply the current values of their estimates. 
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(3.8),3.9) 





Step-1.3 : Compute the new €; by an approximation based on Newton's 


procedure, First, update the value of 


0 
ry = ry - x, B\e" , where Bo =B. (3.10) 


Compute 








Ny; 


C-3 rf 





Feri ie 
iv 3 N, 7 | : 
6 Eijemtr (3.11) 
jel 
; . 
where C=LA,’. 
j=l 
Step-1.4 : Update €; and 8, 
34 = # 


By = §, where k=0,1,....p 


Compute the new value of x, = Day By and iterate the regression, 
starting at Step-J.2, for one more time. 
2. Does the Simple Weibull Regression Adequately Account for the | 
Variability of the Data ? 
The full hierarchical model will be used if the data indicate more 
variability than can be described by the simple Weibull regression. If the data do 


come from the simple Weibull model, then 


U,\e% 
<I j=1,2,.,N, 
By 
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will have thé same distribution as.independent exponential random variables® with 


méan 1. To-assess ‘the variability not accounted for by the:simple regression model, 
“we: compute , 





a2 1 U, ett 
Op = N- 1) reat -1 
WD GT | Hy (3.12) 
M 
and N==N 


Since the variance of a unit exponential is 1, if a, <1 ,then we will stop the model 


fitting with the simple Weibull regression model. If 6% 


x > 1, we go on to the 
hierarchical model to explain the extra-variability. Unfortunately, crude estimate of é, 
uses both the observed and unobserved Uj, A topic for future research is to improve 
the estimate of the variability not accounted for by the simple regression model. 
. Step-1,5 : If 6 $1, stop and use the Weibull regression model only. 
Otherwise, the data will be modeled with the full hierarchical model. A recursive 


procedure to estimate 7 , {B,} and {&} follows and is initialized with {Br} and 
{Es} computed above. 


Si ctl 
f P(U,st}=P\y,(W,)" st}=P(W, s(t Bi) “Vel -e Hi which is the pure Weibull 
distribution. 


ce 





© 
$ 


3.° A Recursive Estimation Procedure for the Parameters of the Full 


Hierarchical Model 





The partial derivative of equation (3.2) with-respect ‘ton, Be , is 
an 


the samé as the derivative of the In-likelihood with respect to n for the model of 


‘Chapter II. The derivative appears in (2.13) and once again the problem is to solve 


the equation f(n)=0. Using the estimates {Bz} and {é;} an initial.estimate for 


n is computed as in equation (2.14) ; that is , 


2 
n° = “fiat (3.13) 


where 


why = -0,5772 


E34 [b,- poe 3.14 
ist) WPv-= 2] (3.14) 
a aVjel 
rr 

DLA 


If nm = <0 ,setn=0.° 

In addition to determining the bounds for a reasonable root to be the 
estimate of n as described in Chapter II, a check for the shape of f(n) is also 
performed. This additional check is performed because two of the simulation 


replications exhibited a very flat f(n). The graph of f(n) for one of the two 


* Same as footnote 4 on page 10. 
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4 ‘, sais 4 = > = ee - ° . 
. 7 4 rad 4. : a ios rs ca ~ ee . ~ oy ae af 7 ° e 
‘¢éplications appears in Figuré-.C 7°. This unusual shape for f(n) tesults in: the 
ne Ro ~ : . 3 





7 


- ‘bisection méthod converging to aitunréasonable root(eg. ff =146 in.theTigure-C). 
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Figure C An Ill-Conditioned Derivative of the In-Likelihood Function for n 


The graph in Figure C indicates that the values of f(n) are asymptotic 
to a fixed nonzero value as n tends to -©; this behavior is the same as in Figure A. 


However, when n is positive, f has very flat curve. The values of f(n) are positive 


0 the graph is one of the only two ill-conditioned cases for f(7) in the simulations. This case ; 
occurred in the gamma/Weibull regression model with O=10, M=15, N=25 and random seed = 
1905665785. The other case occurred with O=15, M=15, N=25 and same random seed. The true value 


of 7 is 1, 














and close-to zér0 but the: curve does not actually ee zero-until 7 is large. The 
following-check is:déne.to detect this situation: If 4 reasonable right-hand bound for 
a ae of (7) cannot be found, ‘check the values.of f(n) and f’(n) for n positive If 
the absolute value of if(n) is ‘less the €, andthe absolute-value of f’(n) is less than 

€3, where: 1 >> €). > €, > O, then the simulation stops for this set of random 
numbers and.new random numbers are generated, Otherwise, perform the bisection 
séarch:as deseribedin. Chapter II. In.the-procedure, £’(n) for n positive is estimated 
as 


[Alnas)~ Flom) 


J(Maw)” ; 


(3.15) 
Tag =°+kxd 


where ns 
Nre-1)™ no+(k-1)xd, 


The algorithm: of this bisection search with the check for flatness of the curve of f 


will be detailed later in Figure D. 


b. The Least Squares Regression Procedure to Salve for {8,} 
The previous estimates of n°, {Br} and {é;} are input to this stage. 


The partial derivative of equation (3.2) with respect to , is 


M 7° 
Onl yy Gre (3.16) 
op, izl op, Ste” op, 


where 
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8 Salrce ae oe) 


2B, jal 














The equation (3.18) can be rewritten as 
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(3.17): 


(3.18) 


(3.19) 


aint eM (=aywe 2 
8, tw <a kal “(B= Bs) Hye ie 


Let: 
- 
and 
A, +w* p A we <. 
c= (Bet Eat MUM yl 
Wy 


The equation (3.20) can be rewritten as 


M Ne Pp 
0=)" fe - Min a 
isl jel he0 


(3.20). 


(3.21) 


(3.22) 


(3.23) 


which are the normal equations for a least squares regression having dependent 


variable z; and independent variables u;,, k=0,1,...,p. 


The solution of (3.23) is 


B=-(UU'WZ 


(3.24) 


where Z and U are defined as in equations (3.25) and (3.26) respectively. 


at 











a ae 





fe “yp 
232 ° « 7 Urry 
‘ Zy ° Urs, 
z=|™ s+ Hany (3.25),(3.26) 
Zul s+ Wap 
Zany ° . ° “ny 





This procedure is an iteratively re-weighted least squares regression; cf. McCullagh 
and Nelder [Ref. 7]. 
c. Newton Procedure to Solve for {&} 


The prior estimates of 7), {8;} ‘ {7 and {é} are input to this 


‘stage. In order to update the values of {€)}, we need to have the partial derivative 


of unconditional In-likelihood function for observer i with respect to €, The 


unconditional In-likelihood function for observer i is the same as equation (2.11). 





‘The derivative is as-follows : 





N; 
ies = >> exp{(y,- x,8°) e*H(-e"*}/y,-x, 8°) : 


i je 


Again, let 


ry = by - Be 


The equation (3.27) can be rewritten as 

















AnL, Ni Cte” 
- —— =h(&) = -C, + ¥r,|-A,+———exp(r,)| . 
0b; jel i Se , 
Ce” 
Treating — as a constant, then 
S,+e" 
PlnL, oh(?; Ce Nj 
= so? a i a ~ Lexp(r,) rj 
0b; i Sire” il 
alnL, ' 
If ee, = h(§) = 0 , the Newton equation for €; is 


alnL, 
a Weil ae He iMfe-Er] - 


* i 


(3.27) 


(3.28) 


(3.29) 


(3.30) 





 - Update the-value of ¢; by solving-the following equation : 





_ go FG) 
pr Si He (3.31) 
ae, 


As noted before, it is possible for equation (3:31) to be badly conditioned. If 
la(é ils £, and ane <e, , where 1>>e,>€,>0, then equation (3.31) is badly 
conditioned. This eee happens occasionally and occurs for only one €; out of 


the M. To avoid this behavior, if lal ils €, and ate <e, , then put 





be aa ; otherwise let €; be the solution to equation (3,31). Return to the bisection 


search procedure and iterate steps 3.a to 3.c until 


“tS 


By Bi Bil 
Br 


i 


ef 


MLTR TR 


Ayu 











< 0.01 for all k=0,..,p and i 

















C. IMPLEMENTATION OF SIMULATION 

Simulation will be used to: study the small sample behavior of the estimators 
of the:five parameters of the gamma/Weibull model. Different numbers of observers 
((M), targets (N,)-and opportunity times (O,) are used as the input to the simulations. 
For each set of input, (O,M,N;), the simulation run starts with the same random seed 
as in Chapter II. The simulation is run until there are 100 (R) estimates of n. A run 
will contain more than 100 replications if the simple Weibull regression model is 
found to describe the variability of the data for some of the replications or if the 
procedure to find the estimate of n is badly conditioned. The statistics of mean bias 
(M.B), mean square error (M.S.E) and their respective standard errors ( S.E(M.B) 
and S.E(M.S.E) ) for all the estimates are computed using the definitions in 
equations (2.16a,b) and (2,17a,b), The averaged uncensoring levels (UC) as 
computed using equation (2.18) are also reported. The statistics for {€;} do not 
include those €; for which |h(é;)|< €,, (€,;=0.1), and |h’(€.)|< €,, (€,=0.01). 

The remainder of this section provides an outline of one replication in the 
simulation, 


1. Input to Simulation 


eM the number of observers ( Possible values are 15 , 25 and 35 ), 


° § the shape parameter of the Weibull distribution for ith observer, 
where i = 1,2,...M; ( Set €; = 0 for alli). 


°7n the parameter for the single-parameter gamma distribution, 
GAM(a,@) , where a = e"; (Setn = 1). 


oN; the number of targets presented to the ith observer , where i = 
1,2,...,.M ; ( Possible values are 15 ,25 and 35). 
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* xX; the k" environmental factor which affects the observation time 
of the i observer whet presented with the j" target. 


¢ ©; — the opportunity time for the i""-observer to détect the j" target; 
(Possible values ate 10-and 15). ie 


e {B,} ‘regression coefficients ae te scale parameter of Weibull 
distribution: , WEI( Mi, a); where k = 0,1,2,. 4p and 


| Do . 
By = cof Ey ;(setp =2; 8 = 08, 8, = -0.2 and 8B, 
ka0 
='0.5 ). 


2. Simulation of Data for the i Observer 
e -Generate a single-parameter gamma r.v..A; from GAM(q,a),having 
density function as in equation (2.1) where a = e" and'n = 1. 


¢ Generate environmental factors rw. x;, ; the x, are independent with 
normal distributions having mean p, and variance o; fork = 
1,2,....p; (Set p = 234, = 1, 0; = 0.5 and p,=2, 0, = 1). : 
Put Xj) = 1 for all ij. 

¢ Compute the Weibull shape parameter p; as in equation (3.1). 


e Generate independent exponential distribution r.v.’s W,, with mean 
1; Le. W,, ~ EXP(1) for all i and j. 


W,,\e! 
¢ Compute the target detection times: U; = Hat : 
i 


° Compute the recorded In-times Y;, , censoring indicator A, and C; 
using the definitions in equation (2. 3) , (2.4) and (2.5) respectively. 
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3. Initial Estimate of £,’s and é; ‘in the Simple Weibull Regression Model 


without the Hierarchical: Gamma r.v.’s 


Step-1.1 : 


Step-1.2 ; 


Step-1.3 : 


Step-1.4 


Step-1.5 





‘Initially, set x8 = nU, and e? =0 


Iteration for regression : 


Compute 1; , wj , uj and z; as in equations (3.3) to (3.6). 


y? y? 


Regress the dependent variable z; against the independent 
variables u;, for {8,} as in equations (3.7) to (3.9). 


Newton’s procedure for estimating €; 


Update the value of 1; as in equation (3.10). 


Compute the new €; as in equation (3.11) with no need to check 
the ill condition. 


Put {é} = &} and {Bj} = {8} for all i=1,2..,.M and 


12 
k=0,1,..,.p. Compute new value of xB = 2 Xi Be and return 
k= 0 


to Step-1.2 for only one iteration. 
Check for extra-variability in the data by computing the sample 
variance 6; for the times Uj as in equation (3.12). 


If a < 1, stop and quote the simple Weibull regression 


model only with the parameters as estimated. 
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4 : 


Otherwise, execute next the recursive procedure. 


4. Recursive Procedure.to Estimate the Parameters of the Full Hierarchical 


Gamma/Weibull Model: 


Input : initial 62 jk=0,1...,p3 &)  ji=1,2,...M. 


Step-2.1 : Initial condition for gamma parameter 7 


Step-2.2 


Pp 
Update x, 8° = © x Bi ; 
ke 0 


Find the initial condition for n using equations (3.13) and (3.14) 


2 
If Tigh, <0 ,setn = 0, : 


Perform the bisection search for n using the algorithm shown in 


Figure D (next page). The algorithm includes a check for the 


flatness of the curve of f(7). 








ect w Ip ability aX. 
Cinrur 


6053 €y0.1 |. 
KEO. 5 €2-0.01} 


J oKeK+1 | 
NEW -Kxé. 
Nee W¥KX8 | 






<TR) * 1(n,) <I> 


(n°) * 1(n,) <9 
Yes 
ne 1° 


nt 
2 


Ni 


: | Teo,)* fn) <> 


Figure D The Algorithm of Bisection Search with Instability Check for 


z =o +kxd 
1 In the Algorithm, we let FI(My)" (aw) sAht where Tray *0 


Naeen=n +(k- 1)xd 


’ 


and k is the current index of increments. 





Iteration for the regression coefficients ,{8,} 


Compute w; , u; and z,; using equations (3.19), (3.21) and 
(3: 22). 


Regress the dependent variable z, against the independent 
variables uj, using equations :(3.24)., (3.25) and (3.26). 


Step-4.1 : Iteration for the In-shape parame tess {é;}, of the Weibull’ 
distribution 


Step-4.2 





. - 
Update x,B° = 2 xp and ry = » y x, Be : 


Compute: H(E;) and aid) using equations (3,28) and 
i 


(3.29) respectively, 
If lH(Er<e, and et alee , for €, = 0.1 and €, = 
i 


0.01, keep &, = é) , Otherwise, update €; using equation 
(3.31). 


: Test for stop criteria 


By Bul [Er 


If on fe 
Bi 


i=1,..,.M, then, stop and record the estimates. Otherwise, set 


M7 NR 
Tu | 











a | < 0.01 for all k=0,..,.p and 








a4 


no=n, By = Pes e = €, and return to Step-2.2. 
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The sitnulation programs for this gamma/Weibull regression model, naméd: 


"SIMULA2", "GAMWEI' and "FVALUE’, are detailed in Appendix D. 


D. RESULTS 

Results of the simulation experiments are presented in Appendix B. The results 
are presented in five sections, Each section shows results for one estimator; the 
gamma parameter n and the Weibull parameters 8, , 8, , 8, and &,. Each section 
contains tables of statistics and histograms of estimate bias for different numbers of 
observers (M or OBS), targets (N or TGT) and opportunity times (O), 

The first two tables in each section present the statistics for different 
opportunity times (O=10 and O=15 in our case ). The next two tables organize the 
information to highlight relations between the number of targets and the number of 
observers and the statistics of the estimates. Histograms of the biases’ of the 

, > estimates are also presented. All parameter estimates appearing in Appendix B are 
from replications for which the full hierarchical gamma/Weibull is estimated. 
Since #1 , the statistics for {€} are for the average of the {€,} for all the 
observers in a replication. The average includes only those €,’s for which the Newton 
procedure converged, The replications for which &, ¢ 1 and the replications in 
which n could not be found are not used to compute the statistics appearing in 
Appendix B; statistics concerning the numbers of these cases appear in Tables A and 


B below. 
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The.statisties in Appendix B use 100 replications. Tables A-and B below show 

, : : the averaged Mngensoring levels, the number of additional replications for which 7 
ae . euis snot be found due to the flatness, Of f, andthe number of times a €,; could not 
be found using Newton procedure due to the smallness of |h(€;)| and |h’(g,)|. The 


results of Tables A and B (next page) are summarized as follows : 


¢ Uncensoring level (UC) : 
A longer opportunity time (O) is reflected in a higher average uncensoring 
level (UC). The opportunity time O = 10 gives about 76% UC, while O 
= 15 gives about 85% UC, This behavior is what we would expect. 


¢ The variability of the data explained by the simple Weibull regression : 
Different parameter estimates in the simple Weibull regression cause the 
different number of cases for which { 6% < 1} for the same input of 
number of observers and number of targets in Tables A and B, In both 
Tables A and B, the cases of { é, < 1} appear more frequently when 
; the number of observers (M) is smaller. This is reasonable behavior since 
more observers tend to provide more evidence that a hierarchical model 
is needed to explain the variability of the data. The manner in which 
6; is computed prevents drawing conclusions concerning the effect of 
censoring on the assessment of unexplained variability. 


e Instability check for 7 : 

There are only two replications in which n could not be found using the 
bisection search for the root of f(n), The function f(n) for one of them 
is plotted as in Figure C, and the function for the other is similar. Both 
replications used the same random seed but had a different set of inputs, 
(M,N;,O;;). These two replications were stopped and new random numbers 
were drawn. This behavior did not occur in the model of Chapter II with 
known parameters for the Weibull distribution. 


e Instability check for € : 
Smaller opportunity times tend to have more replications in which the 
procedure to find €; is badly conditioned. The maximum number for any 
case is 2, 








TABLE A: Comparisons of Uncensoring.levels, Extra Variability Checks and Checks 
of Numerical Instability‘atO =.10 ° _ , 





TABLE B : Comparisons oi’Uncensoring levels, Extra Variability Checks and.Checks: 
of Numerical Instability at O = 15 


Number of} Number of | Average | Number of] Num. unstable 
Observers | Targets uc {{ é, < 1} 


|, (eee 
15 












[35 [ asa [6 [1 


12 All the cases of numerically unstable h and h’ were obtained while replicating the simulations 
starting with random seeds a:999445582, b:444977940, ¢:538790986, d:999445582 & 1534999s0, 
1011786849 and £:444977940. 





In Appendix B, ‘the statistics for the average observer ¢;.do not include those 
individual €, for which the. Newton procedure did not converge. The statistics for A, 
do not inélude those replications for which the simple Weibull regression model was 
used to describe the data. A summaty of the results of the simulation experiment 


shown in Appendix B follows : 


e The estimate of the gamma parameter, 1 

The mean bias tends to be positive. The histograms of biases indicate that _ 
the distribution is skewed to the right. Since the variance of A; ise, the 
positive mean bias of n suggests that for small sample sizes the procedure 
is indicating less between variability than there is in the data. Increasing 
the number of observers tends to decrease the mean bias and mean square 
error. Increasing the opportunity time has little effect on mean bias and 
mean square error. Increasing the number of targets has little effect on 
mean bias and mean square error. As noted in Chapter Il, this behavior 
is reasonable. 


¢ The estimates of the Weibull parameters, 6, through 6, 

The first two tables in parts 2 through 4 of Appendix B indicate that all 
the M.B and M.S.E of {8,} are quite small compared to those for n. This 
is due to the larger sample size used in computing {f,}. An increase in 
the opportun’ty time has little effect on the mean bias and mean square 
error. Increasing the number of targets or observers with a fixed number 
of the other yields small changes in M.B and M.S.E. However, increasing 
both the number of targets and the number of observers does tend to 
decrease the mean bias and mean square error of {8,}. The M.B’s are 
both positive and negative. 


M 
e The average estimate of the Weibull parameter, bre ‘ 
isl 
Since ; =1 in the simulations, the histograms and statistics include the 
averaged &; for all the observers for which the Newton procedure for €; 
converged. Generally, all the magnitudes of M.B and M.S.E of & in the 
first 2 tables of part S in Appendix B are quite small. All the M.B of the 
averaged € are negative. The M.B. and M.S.E. are slightly smaller for . 
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i 


longer opportunity times. Increasing the number of targets or the number 
of observers with the other fixed résults small changes in the M.B and 

MS.E, Increasing both the number .of targets and the number of 
es observers tends to decrease the mean bias and mean square error. 
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IV..CONCLUSIONS AND RECOMMENDATIONS 


The main effort of this thesis is to study the small sample behavior of 
estimators for a hierarchical gamma/Weibull regression model for target detection 
times. This model can also be used-to describe and predict times to failure of similar’ 
machines, (e.g. engines) in different environments. In this thesis, the model assumes 


that. there are M observers. The i'" observer is presented with N, targets, The 


variable U; is the time it takes for observer i to detect the target j. Let Xj , Xj2 ys 


Xjp be the values of environmental variables which. may affect U,. Given: the 


quantified value of the ability for ovserver i, A; = 0, the Uj’s are assumed 
conditionally independent random variables having Weibull distributions, 

WEN, ,€°*) , where the scale parameter has the form 

by = exp{x BI = ox Ea . The random variables {A;} are assumed 
independent from a single parameter gamma distribution, GAM(a,a) witha = e". A 
numerical procedure based on maximum likelihood is used to estimate the 
parameters of the model. The numerical procedure is iterative and uses a bisection 
root search method for estimating the gamma parameter n , least squares regression 
to estimate the Weibull In-scale parameters 8) through 8, and an approximate 
Newton procedure to estimate the Weibull In-shape parameters {¢,}. Simulation is 
used to study the behavior of the estimates for small sample sizes. Generally, the 


numerical estimation procedures work well. There are few numerical problems. All 
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‘the mnean biases.and mean square errors of the estimates tend to be small compared: 
to the theoretical inputs-to the simulation.. 

Topics for future research include replacing the one parameter gamma 
distribution with the two parameter gamma distribution in the hierarchical model and 
exploring the use of EM methodology to estimate model parameters. The two 
parameter gamma distribution arises as the posterior distribution of the current 
model. Another area is to study using the fitted hierarchical gamma/Weibull model 
to predict future performance. 

It is hoped that the hierarchical gamma/Weibull regression model will be a 
useful tool to describe and predict one aspect of the effect of human: performance 


on the battlefield. 


43 





‘LIST OF REFERENCES 


Gaver, D, P., Jacobs, P..A.,.and.O’Muircheartaigh, I., Hierarchical Models for 
Target Detection: Times, Naval Postgraduate School Technical Report. To 
appear. 


Bae, Deok Hwan, Models for Target Detection Times, Naval Postgraduate School 
Thesis, September 1989. 


Johnson,L. W. and:Riess, R. D., Numerical Analysis, Addison-Wesley,Inc. New 
York, 1977, 


Peter A, Stark,Introduction to Numerical methods, Macmillan Publishing Co.,Inc. 
New York ,1970, 


Lewis, P, A. and Uribe,L.,7he New Naval Postgraduate School Random Number 
Package LLRANDOM II, Naval Postgraduate School Technical Report, NPS55- 
81-005, Monterey, CA, 1981. 


Leonard Gilman and Allen J. Rose, APL An Interactive Approach, Third Edition, 
JOHN WILEY & SONS,Inc. New York, 1984. 


P, McCullagh and J.A. Nelder, Generalized Linear Models, Chapman and Hall, 
New York, 1983. 


44 








APPENDIX A. STATISTICS AND GRAPHICS FOR THE ESTIMATE OF 
4, IN THE SINGLE PARAMETER GAMMA MODEL 


“Table 1. Mean Bias, Mean Square Error, Std. Error and Uncensored Level for 7 at 
O = 10 inthe Gamma Model 


Number of Number of 


el Sere 
| Observers | Targets [ MBGE) | _MSEGE) | 
is saat nisi os — i 
nF aH 
0.1003 aa 5905 


0.0566 (0.0304) 0.0935 (0 012 ) 


0.0682 (0.0274) 0.0784 (0.0108) 


0.058 (0.0281) 0.0807 (0.0129) | 82° 
28 0.0088 (0.0253) | 0.0626 (0.0092) 


0,05 (0.0103) 





Table 2, Mean Bias, Mean Square Error, Std. Error and Uncensored Level for 7 at 
O = 15 in the Gamma Model 


Number of | Number of 
Observers Targets 
0.1886 (0.0413) K9% 
0.1875 (0.026 ) 
0.053. (0.0452) 0.2032 (0.052 ) 
0.094 (0.0505) 
0.0517 (0.03 ) 0.0911 (0.0115) 90% 
0.0706 (0.0274) 
0.0519 (0.0277) 0.0781 (0.0132) x 
0.0085 (0.0246) 0.0593 (0.0086) 


0.0471 (0.0216) 0.0481 (0.009 ) 








Table 3. Tétidencies of Mean Bias and Mean Square Error for y at O = 16.in the 
ne Gamma Model. 






M.B(MS.E) of #7 











{1s oBs_| * 0.073(6.248)_ | __0.1030.197)_ | 0.067(0.217)__ | 
| 25-0BS | 0-1. (0.105) | 0.057(0.094) | 0.068(0.078) | 
{35 0Bs | 0.058(0.081)__ | __0.009(0.063)__|__0.087(0.05) | 








‘O = 10 ( UC ; 82 PERCENT ) 
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RAS OF ETA 


N=15 N=25 N=35:M=15 M=25 M=35 


Figure 1, The Tendencies of the Bias of y with changes in OBS and TGT at O= 10 
in the Gamma Model: M = number of observers ; N = number of 
targets. 
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Figure 2. The Tendencies of the Bias of 7 with changes in OBS and TGT at O= 10 
inthe Gamma Model: M. = number of-observers ; N = number of 
targets. 
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Figure 3. The Tendencies of the Bias of 7 with changes in OBS and TGT at O= 10 
in the Gamma Mode: MM = number of observers ; No = number of 


targets. 
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Figure 4. Comparison of the Bias of 7 between Different numbers of Observers with 
TGT = 15 at O=10 in the Gamma Modek TGT = the number of 
targets. 
rd 
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Figure 5, Comparison of the Bias of between Different numbers of Observers with 
TGT = 25 at O=10 in the Gamma Modek TGT = the number of 


targets, 
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Figure 6. Comparison of the Bias of 7 between Different numbers of Observers with 
TGT = 35 at O= 10 in the Gamma Model: TGT = the number of 
targets. 
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Table. 4, Tendencies. of Mean Bias and Mean Square: Error fory at O = 15 in the 
sa Model 






M: BM. S.E) of # 










0:061(0.189), .085(0.158) 0.053(0.203) | 
25 OBS 0.094(0.1 ) ~0.052(0.091) 0.071(0.078) 
| 35:OBS 0.052(0.078)_ 0.008(0.059) 0.047(0.048) 














O= 165 ( UC: 90 PERCENT } 
M= 15 <¢=> N= 15 


EWS GF ETA 






n2anater? pentahbanaad ee ceree 







2 4 & 
N=15 N= 25 N= 35:M=15 M=25 M= 355 


Figure 7, The Tendencies of the Bias of 7 with changes in OBS and TGT at O= 15 
in the Gamma Model: M = number of observers ; N = number of 
targets. 
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Figure 8. The Tendencies of the Bias of y with changes in OBS and TGT at O= 15 
in the Gamma Model: Mf = number of observers ; N = number of 
targets. 
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Figure 9. The Tendencies of the Bias of 7 with changes in OBS and TGT at O= 15 
in the Gamma Model: M = number of observers ; N = number of 
targets. 
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Figure 10. Comparison of the Bias of 7 between Different numbers of Observers 
with TGT = [5 at O= 15 in the Gamma Model; TOT = the number 
of targets. 
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Comparison of the Bias of y between Different numbers of Observers 
With TGT = 25 at O=.15 in the Gamma Model: TGT = the number 
of targets. 


Figure {1 
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Figure 12, Comparison of the Bias of 7 between Different numbers of Observers 
with TGT = 35 at O= 15 in the Gamma Model: TGT = the number 
of targets. 
4 


35 








APPENDIX B. STATISTICS AND GRAPHICS FOR ALL THE 
ESTIMATES IN GAMMA/WEIBULL REGRESSION MODEL 


1. eee FOR y 


Table 5.. Meat Bias, Mean Square Error.and Std; Error for 7 at oO = 10 in the 
GAM/WEI Regression Model 


| Number of Number of 
Observers | Targets [| MBSE) | MS.ES.E) 
~_3299(,0892) 
pe 8573(.0455) | 2274(,0445) | 
1074(.0352 .1328(.0216) 


0658(.0307) 0968(.0129) 


0577(,0070) 






























Table 6, Mean Bias, Mean Square Error and Std. Error for 7 at O = 15 in the 
GAM/WETI Regression Model 



















Number of | Number of 















Observers Targets M.B(S.E) M.S.E(S.E) 
| .2107(.0561) 3524(.1023) 
25 .1934(.0451) .2365(.0458) 








35 1117(.0442) 2038(.0319) 
1170(.04 18) 1851(.0454) 


is | 18500356) 1315(.0238) 
35 0982(.0133 
_9566(,0072) 
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Table 7. Tendencies of Mean Bias and Mean Square Error for 7 at O = 10 in tlie 
GAM. iW El Regression. Model 


. fo . ue M. B(M.S:E) of 
i 35 

‘ roa sya 076(.196) 
“OHS 057(.100) .142(.123) 


.066(,097) .094(,058) 
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Figure 13. Comparison of the Bias of 7 between Different numbers of Observers 
with TGT=15 at O= 10 in the GAM/WEI Regression Mode TGT 
* = the number of targets. 
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Comparison of the Bias of y between Different numbers of Observers 
with TGT=25 at O= 10 in the GAM/WEI Regression Model: TGT 
= the number of targets. 
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Figure 15. Comparison of the Bias of 7 between Different numbers of Observers 
with TGT= 35 at O= 10 in the GAM/WEI Regression Model; TOT 
= the number of targets. 





Table 8. Tendencies of Mean Bias and Mean:Square Error for 7 at O = 15 in-the 
GAM/WEI Regression Model 


a _ M.B(MS.E) of # ne | 
1 fs TCT 25 TGT 35TGT 
fl _211(,352) 
“117(.185) 
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Figure 16. Comparison of the Bias of » between Different numbers of Observers 
with TGT= 15 at O=15 in the GAM/WEI Regression Model: TGT 
= the number of targets. . 
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Figure ge Comparison ‘of the Bias of 4 between Different numbers of Observers 
with TGT=25 at O= 15 in the GAM/WEI Regression Model: TOT 
= the number of targets. 
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TGT = 35,0 = 15 ( UC: 85: PERCENT ) 
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Figure 18. Comparison of the Bias of 7 between Different numbers of Observers 
With TGT=35 at O= 15 in the GAM/WEI Regression Model: TGT 
= the number of targets. 





2, TABLES AND HISTOGRAMS FOR 0 


Table 9, Mean Bias.Mean Square Error and Std. Error for £0 at O = 10 in the 
GAM/WEI Repression Model a 









af Nuinber of | 


; ~ Number of. 
a Observers 


Targets 


ps ; aN, 0317) TON 0123) 
: 0034 0278) 9764 0109) 










a “OIG. 0222) 0487 0032) 


+,009-4(.0197) .0383(.0056) 


-.0067(.0198) 0386(.0067) _ 
-.0398(.0176) 0319(.0047) 
| 0062(.0152) 0227(.0027) 


Table 10, Mean Bias,Mean Square Error and Std, Error for fQ at O = 15 in the 
GAM/WE1] Regression Model 





= 

















Number of 

Targets M.B(S.E) 
~.0609(,0316) 
-.0234(.0271) 
.0263(.0224) 
.0365{.0244) 
-.Y292(.0210) 


Number of 
Observers M.S.E(S.E) 
1015(.0127) 
.0725(.0095) 
.0498(.0063) 
.0595(.0073) 
V439(.0049) 




































335 -.0185(.0184) .0335(.0033) 
| 1S -.0043(.0191) .0358(.0066) 
25 -.0454(.0176) .0324(.00-41) 





-.0021(.0138) .0187(.0023) 








“ 


Table I1.. Tendencies-of Mean Bias and. Mean Square Error for {0 at O = 10 inthe 
: a GAME] Regression Model: 


MOISE) of hy 


23 OBS 46(.051) -.009(.038) 
rBa00032) 006(.023) 


= 10 (UC: 76 aha ) 
= 25 
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Figure 19, Comparison of the Bias of £0 between Different numbers of targets with 
OBS = 15 at O= 10 in the GAM/WEI Regression Model: OBS = the 
number of observers 
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' Figure 20. Comparison of the Bias of f0.between Different numbers of targets with 
OBS = 25 at O= 10 in the GAM/WEI Regression Model: OBS = the 
number of observers 
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Figure 21, Comparison of the Bias of £0 between Different numbers of targets with 
BS = 35 at O= 10 in the GAM/WEI Regression Model: OBS = the 
number of observers 





Table 12. Tendencies-of Mean Bias and Mean Square Error for £0 at'O = 15'in the 
_ GAM/WEI Regression Model. : ; 

-_ M.B(MS.E)-of f, . 
Hie ster ster] 
{15 ops” | -.061(.102)_ | =.023(,073)_ | .026(.050) | 
[25.0Bs. | .037(.060) | -.029(,044)_ | -.018(.033) 
L3soBs_ | -.004(.036) | -.045(.032) -.002(.019) 
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Figure 22. Comparison of the Bias of {0 between Different numbers of targets with 
OBS = 15 at O= 15 in the GAM/WEI Regression Model: OBS = the 
number of observers 
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Figure 23, Comparison of the Bias of £0 between Different numbers of targets with 
OBS = 25 at O= 15 in the GAM/WEI Regression Model: OBS = the 
number of observers 
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Comparison of the Bias of £0 between Different numbers of targets with 
OBS = 35 at O= 15 in the GAM/WEI Regression Model: OBS = the 


number of chservers 


Figure 24, 
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3.. TABLES.AND: HISTOGRAMS FOR f1 


Table 13. Mean Bias;Mean Square Error and Std; Frror for B1 at O = 10 in the 
GAM/WEI ‘Regression Model 


_ Number of i: See 
“Observers | "Targets 
{is _[__o1soors) | ___.0243(.0032)__| 
Is | 5 | _-0002(.0140) | 0193/,0026)_ 
1S | =.0154(.0132) | .0173(.0026) 
25 | 5 -.0006(:0096) | 091.0014) 
p35 -.0036(.0080) | _-0063(.0010)_ | 

35 ACH SIGHT 












Table 14. Mean Bias,Mean Square Error and Std, Error for #1 at O = 15 in the 
GAM/WEI Regression Model 


Number of | Number of 


Observers | Targets Mass a 
pS | 007960145) | 0205(.0028) 

35 orca ame 0014) 

| is -.0214g.oi27) | _0162(.0024) 
25 | 25 | _.0036(,0090) | 009.0013) 
0086, 0073) 
25 
35 000940067) 





Table 15. Tendencies.of Mean Bias and Mean.Square Error for f1 at O = 10 in the 
; GAM/NEL Regression Model . ; 

M. au S.E) of Be 
: 













OBS = 15, O = 10 ( UC: 76 PERCENT ) 
N= 25 


C 
BUS OF BETAY 


0 
BUS OF BETAL 





Figure 25. Comparison of the Bias of £1 between Different numbers of targets with 
OBS = 15 at O= 10 in the GAM/WEI Regression Model OBS = the 
P number of observers 





OBS = 25:, 0 = 10-( UC : 76 PERCENT ) 


N. = 25 





Figure 26. Comparison of the Bias of £1 between Different numbers of targets with 
OBS = 25 at O= 10 in the GAM/WEI Regression Model: OBS = the 
number of observers 





76 PERCENT ) 
N-= 25 





Figure 27. Comparison of the Bias of #1 between Different numbers of targets with 
OBS = 35 at O= 10 in the GAM/WEI Regression Model: OBS = the 
number of observers 





Table 16. Tendencies of Mean Bias.and Mean a Error for B1 at O = 15 in the 
oe GAM/WEI Regression. Model 














M. BM S. E) of B, 


| 
CET eC 
[25 ons] -02i(016) | 004.008) | 0054.06) 
[35 ons 10150013) eee | ,001(,004) 









= 15,0=15( UC: 85 PERCENT ) 
N = 25 


8 


BAS OF BETAY 


8 


0 O58 
BIAS OF BETAL 





’ Figure 28. Comparisen of the Bias of £1 between Different numbers of targets with 
OBS = 15 at O= 15 in the GAM/WEI Regression Model: OBS = the 
number of observers 








OBS"= 25 , 0: = 15 ( UC.: 85 PERCENT ) 
N= 15 


_* 


[) 
DUS OF BETA! 





Figure 29, Comparison of the Bias of £1 between Different numbers of targets with 
OBS = 25 at O= 15 in the GAM/WEI Regression Model: OBS = the 
number of observers 








ok OBS = 35 , 0 = 15 ( UC : 85 PERCENT ) 
; ; A NS 15 ‘N =.25 





Figure 30. Comparison of the Bias of #1 between Different numbers of targets with 
OBS = 35 at O= 15 in the GAM/WEI Regression Model: OBS = the 
number of observers 








‘4, FABLES AND. HISTOGRAMS FOR.2 
Table 17. Mean. Bias;Mean Square Error and Std. Error-for £2 at O = 10 in the 
Par “GAM/WEI Regression Model —_ : 

' Number of 
| Observers. 



















_areets [MBG E) | MS.ES.E) | 
[15 | .0167(,0097)_ | _.0096(:0012) | 
{3s | -.0097(.0053)_ | .0029(.0004) 
2 |__| .0040(,0054) | 0029.00) 


Table 18. Mean Bias,Mean Square Error and Std. Error for 82 at O = 15 in the 
GAM/WEI Regression Model 


Number of | Number of 
ial 
.0185(,0091) 
.0055(,0070) 
-.0124(,0052) 
~-0015(,0069) 
.0072(.0050) 
_-0024/,0042) 

— _0016(.0057) 

35 2 0086(.0045) 
-.0007(.0032) 


















































Cry 


wn 





Table 19, Tendencies of Mean Bias and Mean Square Error for £2 at O = 10 in-the 
GAM/WEI Regression Model 

2 M.B(MAS.E) of B; | 
7 | stor | 2s TGT | 5 TCT 
{_25.0BS_|_ -.005(.005)_ | __.004(.003)_|—.002/.002) 











OBS = 15, 0 = 10 ( UC: 76 PERCENT ) 
N= 15 N = 25 


2 9 02 
BUS OF BETAZ 





Figure 31. Comparison of the Bias of £2 between Different numbers of targets with 
OBS = 15 at O= 10 in the GAM/WEI Regression Model: OBS = the 
number of observers 





OBS = 25, 0-= 10 ( UC: 76 PERCENT ) 





Figure 32. Comparison of the Bias of £2 between Different numbers of targets with 
OBS = 25 at O= 10 in the GAM/WEI Regression Model: OBS = the 
number of observers 





OBS = 35 ,.0 = 10 ( UC : 76 PERCENT ). 
N 


N-= 15. 2°25 





Figure 33. Comparison of the Bias of £2 between Different numbers of targets with 
OBS = 35 at O= 10 in the GAM/WET Regression Model: OBS = the 
number of observers 
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Table 20. Tendencies of Mean Bias and Mean Square Error for f2 at O = 15 in the 


‘ore GAM/WEI. ou Model 


a 
{_250ns | .001(005) | 007.002) [002 002) | 














,O = 15 ( UC : 85 PERCENT ) 
N = 25 


BIAS OF BETAZ 





Figure 34. Comparison of the Bias of £2 between Different numbers of targets with 
OBS = 15 at O= 15 in the GAM/WEI Regression Model: OBS = the 
number of observers 





OBS =°25, 0 = 15 (UC: 85 PERCENT ) 
Nats : N= 25 


0.2 ° 0 G2 
BIAS OF BETA2 





Figure 35. Comparison of the Bias of $2 between Different numbers of targets with 
OBS = 25 at O= 15 in the GAM/WEI Regression Model: OBS = the 
number of observers 
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Figure 36, 








OBS. = 35 , 0.= 45 ( UC: 85 PERCENT ) 
N = 25 


Comparison of the Bias of £2 between Different numbers of targets with 
OBS = 35 at O= 15 in the GAM/WEI Regression Model: OBS = the 
number of observers 











5. TABLES AND-HISTOGRAMS FOR & 


Table 21. Mean Bias, Mean Square Error and Std. Error for € ‘at O = 10 in the 
GAM/WEL Regression Model 









| Number of | Number of sig 
See 
. 
RET 
1s | 0226¢.0058) [00830008 
25 0013(.0002) 


: 0025. 0038) .0014(,0002) 
ne: 


Table 22, Mean Bias, Mean Square Error and Std. Error for § at O = 15 in the 
GAM/WEI Regression Model 


Number of 
Targets M.B(S.F) 
—— 
Pts | -.o1siy.o0a7) | 0025(,0003) 
25 
Serre -.0192(.0041) 
<n ee ee ~-0049(,0033) 0011(.0002) 
| 35 | -.0083(.0025) | _—-0007(.0001) 

























-.0196(.0044) 0023(.0003) 
-.0065(.0036) 0013(.0002) 


-.0071(.0028) 0008(.0001) 






Number of 
Observers 
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Table 23. Tendencies of Mean Bias-and Mean Square Error for € at O = 10 in the 
_ GAM/WEI Regression Model 


a 2 e M.B(M.S.E) of 
-  - 45 TGT 25 TGT 35 TGT 
; 15 OBS -.032(.005) -016(.003) | —_-.013(.002)__| 


TCT = 15, 0 = 10 ( UC: 76 PERCENT ) 
M = 15 M = 25 


“OA = 


2 i} 
BUS OF x! BIAS OF XI 





Figure 37, Comparison of the Bias of ¢ between Different numbers of Obse rvers 
with TGT=15 at O= 10 in the GAM/WE! Regression Model: TGT 
‘ = the number of targets. 
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TGT = 25,0 = 1d ( UC : 76 PERCENT ) 
M = 15 ‘M = 25 





Figure 38. Comparison of the Bias of ¢ between Different numbers of Observers 
with TGT=25 at O= 10 in the GAM/WEI Regression Model: TGT 
= the number of targets. 
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Tet = 35, O = 10 ( UC: 76 PERCENT ) 
M = 15 M = 25 
8 


Figure 39. Comparison of the Bias of ¢ between Different numbers of Obse rvers 
with TGT= 35 at O= 19 in the GAM/WEI Regression Model: TGT 
= the number of targets. 














Table 24. Tendencies of Mean-Bias and Mean ei Error for € at O = 15 in the 
_GAM/WEI Regression Model 


) 
re ee 
[35 ons [019.003 | -005¢008) [-008¢001) 














TGT = 15, 0 = 15 ( UC: 85 PERCENT ) 
M = 15 M = 25 





Figure 40. Comparison of the Bias of € between Different numbers of Observers 
with TGT=15 at O=15 in the GAM/WEI Regression Modek TGT 
= the number of targets. 





Figure 41. 


TCT = 25, 0 = 15-(.UC : 85. PERCENT ) 





Comparison of the Bias of ¢ between Different numbers of Observers 
with TGT=25 at O=15 in the GAM/WEI Regression Model: TGT 
= the number of targets. 
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PERCENT ): 
M = 25 












Comparison of the Bias of ¢ between Different numbers of Observers 
with TGT=35 at O= 15 in the GAM/WEI Regression Model; TGT 
= the number of targets. 


Figure 42. 


APPENDIX C. SIMULATION PROGRAM FOR THE ESTIMATE OF 4 


3 RRR 
[1] 

[2] 

[3] 

C4] 

[5] 

[6] 

[7] 

C8] 

[9] 

[10] 
[11] 
[12] 
[13] 
C14] 
C15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
C24) 
[25] 
[26] 
[27] 
(28) 
[29] 
[30] 
[31] 
[32] 
[33] 
C34] 
[35] 
[36] 
(37] 
[38] 
[39] 
£40] 
C44] 
[42] 


IN THE SINGLE PARAMETER GAMMA MODEL 


V ReSIMULAI3M:N3sO3I3sd3K3L3ETA3OUT;PAGE;AVUCL3MB;SEMB;MSE;SEMNSE 


nARMAAARAAA KEY GLOBAL VARIABES : BOOK1 ; BOOK2 ; UCL mnnnannnaa 
AAAARAAARAAAAAAARAAAAAAAAAAR DATA INPUT AARAAAAARAARAARAAAARAAAAA 
M<« 15 25 35 

N+ 15 25 35 

0+ 10 15 

ETA¢1 

L<J<K<L<0 

AARARAARAAAAAAARAAA START LOOPING THE SIMULATION 1 naaRARARAAAAAA 
LOOPI:I[<I+1 
LOOPS 3Jed +1 
LOOPK: K<K+1 

AAARARAARAAAARARAAAR RESET THE RANDOM NUMBER SEED aAnAAAAARAARA 
O«OCI) ,ML7],NCKI 

ORL<+466801743 
LOOPL: L<L+1 

BARARAAAARAAAARAA CALL GAMMA ESTIMATE FUNCTION RAARAAAARAAAAAARAA 
RRR<ORL 

OUT<+GANMA (OCI],MC7I,NCK]) 

AAAAAAAR RECORD UCL AND ALL THE ESTIMATES IN REPLICATIONS arnanaa 
*((L=1),L>1)/INTZ1,G01 
INTZ1:PAGE<(1,pOUT)p (OUT-0,ETA) 

>+LOOPL 

GO1:PAGE<PAGE,(1] (OUT-0,ETA) 
TESTL: 

+(L<100)/LOOPL 

ARARAAR COMPUTE AVE. UNCENSORED LEVEL AND ALL THE STATISTICS an 
AVUCL<( (+/PAGEL31] +L) 

MB<(+/PAGECL32])¢L 

SENB+((+/ (PAGEL 32] -MB)*2)#(Lx2-2))x0.5 

MSE<(+/PAGE(L32]*2)4+L 

SENSE<((+/( (PAGEL 32] *2)-MSE)x*2)4#(LxL-2))*0,5 

L<0 

AAAAAARAAR KEEP ALL THE INFORMATIONS FOR GRAPHICS AND TABLES ana 
>(((Ix7xK)=1), (2x7xK 41) /INTZ2 ,GO2 

INTZ2: 

BOOK1i«(1,pPAGEL 32] )pPAGEL; 2] 

BOOK2<(1,4)p (MB,SENB,NSE,SEMSE) 

UCL+AVUCL 

>TESTIK 

G02: 

BOOK1+B00K1,(1] PAGEC3;2] 

BOOK2+BOOK2 , (11 (MB,SEMB MSE ,SEMSE) 


 [ysy 


UCL<UCL,AVUCL, 


E4u) TESTK: 


Cus. 
C46]: 
[a7] 

- [4a] 


[49] 


>(K<pN)/LOOPK 


Keo 
+(J<pM)/LOOPS 
J*0: 


>+(I<p0)/LOOPL 
V 





[42 
[2] 
[3] 
C4] 
[5] 
16] 
[7] 
[8] 
[9] 
[10] 
[11] 
[12] 
[13] 
C14] 
[15] 


C16]. 


[17] 
C18] 
[19] 
[20] 
(C21) 
[22] 
[23] 
t24) 
[25] 
C26] 
(27) 
[28] 
[29] 
[30] 
[34] 
[32] 
[33] 
C34] 
[35] 
[36] 
[37] 
[38] 
[39] 
C40] 
C41) 
C42] 
C43] 
C44) 
C45] 
C46) 
C47] 
[4s] 


V R*GANMA INPUT3M3N;3O03ETA$ETAO$XI3XIO3ASWW3MU3U3Y3DELTAsC3XB3RR 
3M15M2;7;DATA;K;FL3FM3FR;LB3;MP3RB3S 

AAAARRAAAAAA SIMULATION INPUTS AND THEORETICAL PARAMETERS aanan 
M+«INPUT (2). 

N+INPUT(3] 

O+(M,N)pINPUTE1] 

ETA+1 

MU<u, 2 
XI0¢Q(N,M)pXI<Mp0 

RmARARAAAAAAAAAAAAR RANDOM NUMBER GENERATIONS ananaarnananaranrna 
A<&(N,M)pM GAMRAND((xETA),+*ETA) 

WW<(M,N)p (MXN). EXPRAND 1 

RAAARAAAAAAAR DETECTION TIME AND UNCENSORING DATA aaRAARAAAAAAA 
U+MUx (WHA )x (xXT0) 

A«WW+10 

Y+(@U}L (@0) 

C+t+/[2] DELTA+(eU )s(@0) 

0+10 

XB¢OMU 

AAAAAAAAAARA INITIAL CONDITION FOR GAMMA PARAMETER ,ETA ananaan 
RR (Y-XB)x*(71xXI0) 

M1«70.5772 

M2¢(+/+/DELTAXRR*2 )t.(+/C) 

ETA0<0 

+ (0>Te(( (01 x2 #76 )+(M2%2 )~(M1%2))/BISECTION 

ETAO+” 1x0 

ARAAA FIND BOUNDS AND BISECTION SEARCH FOR GAMMA PARAMETER ,ETA aa 
BISECTION: K+0 

DATA (S<+4+/[2]*RR),C 

FN<ETAO FVALUE DATA 

AAAARAARARBAAAR FIND THE RIGHT AND LEFT BOUND FOR ETA aananananaan 
BOUND: K*K+1 

FL+(LB+ETAO-Kx0.5) FVALUE DATA 

FR<(RB«ETAO+KK0.5) FVALUE DATA 

+ (((xFM)x (xFL) )<0)/SELRB 

>(((xFM)x (xFR) )<0)/SETLB 

+BOUND 
SELRB: RB«ETAO 

+BISECT 

SETLB: LB¢ETAO 

AARRAAARAAARAAAAR PERFORM THE BISECTION SEARCH FOR ETA ananannanan 
BISECT: 

FM+«(MP<(LB+RB)+2) FVALUE DATA 

FL+LB FVALUE DATA 

>((RB-LB)<(|0,01xMP))/END 

>(CC(XFM)x (XFL) )<0), (((xFM)x (xPL))>0))/SWAPR,SWAPL 
SWAPR:RB<MP : 

+BISECTL 

SWAPL: LEMP 

>*BISEC? 
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[C49] maaRARAAAAAAAAAAARARAAA OUTPUT THE ESTIMATES aanaaaranaAAAnaAAA 
> 7° 50) END: 2 
[51] ETAO<(RB+LB)+2 | 
[52] .R<((+/C)+#MxN) ,ETAO. 
Vv 


ae 


V RETA FVALUE DATA3M3S;3C3SUBF;I a 
C1] Me(pDATA)+#2 
[2]  S«MADATA 
[3] C<MVDATA 
[4]  SUBF+0 
[5] I<0 
C6] LOOP:I<«I+1 
C7) >(CC(I]=0)/CHECKI 
[8]  SUBF«SUBF+(+/+(xETA)+(1CLI])-1) 
[9] CHECKI: 
[10] +(I<M)/LOOP 
C14] ReCvETA)xSUBF+(+/ETA+((S-C)#S+xETA)+°1x@(S+*xETA)) 
V 








APPENDIX D. SIMULATION PROGRAM FOR THE ESTIMATES OF 
ALL. PARAMETERS.IN GAMMA/WEIBULL REGRESSION MODEL 


3; SEMB-;MSE; SENSE; RRR 


C3] M+ 15 25 35 
C4] Ne 15 25 35 
£57 O+ 10 15 


C6] ETA«1 
[7] B0+0,.8 
[8] B1«"0,2 
£9] B2«0.5 
[10] XI+0 


(12}) T«d+K+L<E+X+0 
[13] LOOPI:I+I+1 
C14) LOOPI:J+J+1 
(15) LOOPK:K«K+1 


* £18] ORD<+466801743 

(19) LOOPL:L+L+1 

£20] apmrmrmAAnAAAAAR CALL GAMMA/WEIBULL REGRESSION FUNCTION 
* [21] RRAR<DRL 


C24] +(OUT[2]*~99999)/CHECKETA 
(25] E+EH+1 

[26] L<L-1 

{27] +ZESTL 


[29] CHECKETA: 

[30] +(OUTC2]4#~100000)/G0GO 
(31) Q«'BAD ETA! 

[32] O+«RRR 

C33] LeL-1 

C34] +TESTL 


[36] GoGo: 
[37] +((5=1),L>1)/IN72Z1,C01 
- [38] INTZ1:PAGE<(1, (pOUT)-1)p(("140UT)-0,ETA,BO,B1,B2,XI) 
[39] X<X+("1+0U7) 
C40] +ZO0OPL 
C41] GO1:PAGE<PAGE,[1]((°1¥0UP)-0,ETA,BO,B1,B2,XI) 
[42] X<«X+(71+0UT) 






























V ReSIMULA2 MsNsOsIsJsKsLsXsEsETA; ;BO$3B1:3B2;XI;0UT; PACE; AVUCEX 3MB 
C43) mamnarnaaannan KEY GLOBAL VARIABES : BOOK1 ; BOOK2 :; UCLEX arannannaa 


[2] ARAAARARRAAAARRARAARAAAAAAA DATA INPUT RAR ARAARRAAARAARAARAAAA 


[14] marnaarnAAAAARARAR START LOOPING THE SIMULATION araaanAARARAARAAR 


[16] pAAAAAAARARAARARARA RESET THE RANDOM NUMBER SEED nanannaarRaAanananana 


AARRAARARAR 


[23] naananananaaa SKIP THE INSUFFICIENT VARIABILITY CASE FOR GAMMA 


[28] anarnaaR RECORD THE RANDOM SEED FOR THE ILL CONDITION OF ETA 


C35] ananARAARA RECORD UCL AND ALL THE ESTIMATES IN REPLICATIONS 











| ARAARAAA RECORD THE RANDOM SEED FOR THE ILL CASE OF XI parannann 


+(.(7140UT )=0) /TESTL: 
O<+'BAD XI! 
O«RRR 


. APARAARARARRAARAARARAPARAAAARAARAAARARARAAAAAAARARARARAARAANARARAARRA 


TESIL: . 
+(Z<100)/LOOPL 


—armarnaaR COMPUTE AVE. UNCENSOFP* ) LEVEL AND ALL THE STATISTICS a 


AVUCEX+( (+/PAGE( 31] )#Z),E,X 

MB+(+/[1] PAGE*+ 0 1 ¥PAGE)+L 

SEMB+ ((+/[1] (PAGE~ (pPAGE )pMB)*2)4(LxL-2))*0.5 
MSE«(+/(1] PAGE*2)4+L . 

SEMSE+( (+/01] ( (PAGE*2 )- (pPAGE) pMSE )*2)4(LxL-2))*0.5 
LE+X+0 


* mmaaannAA KEEP ALL THE INFORMATIONS FOR GRAPHICS AND TABLES aaa 


+(((IxIxK)=1), (Ix7xK)#1)/INTZ2 ,GO2 
INT22: 

BOOK1«(1,pPAGE)pPAGE . 
BOOK2+(1,4,pMB)p(NB,SEMB ,MSE,SEMSE) 
UCLEX« 1 3 pAVUCEX 

+TESTK 

G02: 

BOOK1«BO0OK1,(1] PAGE 

BOOK2«BOOK2 ,(11]((4,pMB)pMB,SEMB,MSE,SEMSE) 
UCLEX<UCLEX,({1] AVUCEX 

TESTER: 

+(K<pN)/LOOPK 

K+0 


>(J<pM)/LOOPI 
J<0 
+(I<p0)/LOOPI 
v 








VGAMWEICTO]V . 

V ReGAMWEI INPUT3M3N30;ETA3ETA03B3BO03B13B2sXIs3XI0OsAsWW3X13X2 
3MU3U3Y3DELTA3C3XB3RR3W3U03U130U23UU3Z3T3VAR3M13M2;7;DATA3K 
:FL;FL03FM;FR3FRO;LB;MP3RB3CS3H; DE; DIFFXI;BK3S3J3X$XX3BADXI 

f [1] anarananaaaa SIMULATION INPUTS AND THEORETICAL PARAMETERS anannana 
C2] M+INPUTLE2) 
[3] N<INPUTL3) 

8 C4] 0«(M,N)pINPUT(1] 


C5] ETA<1 
[6] BO0«0.8 
£7]. B1+70.,2 
£8] B2«0.5 


(9) XI0¢R (NM) pXI<Mpo 
[10] aAaanaArmanAAAAARAA RANDOM NUMBER GENERATIONS anarnanAAAAAARAAAAAAA 
£11) A+Q(N,M)pM GAMRAND ((xETA),+xETA) 
£123) WW<(N,N)p (MXN) EXPRAND 1 
£13] Xa<(N,N)p (MXN) NORRAND 1 0.5 
C14] .X2«(M,N)p(MXN) NORRAND 2 1 
£153 AAAAARARAAAA DETECTION TIME AND UNCENSORING DATA nanAAARARAAARAAR 
C16] MU+* (BO+(B1xX1)+B2xX2) 
[417] U<MUX (WH#A)* (xXI0) 
C18] AWiW<+NU+10 
£19] Y<(ey)l (eo) 
(20] . C#+/£2] DELTA«(@U)s$(@0) 


[21] 0<10 

€ [22] AAAAAAAAARA INITIAL ESTIMATES WITHOUT GAMMA R.V. INVOLVED naanaaan 
[23] XB<oU 
C24] I<0 


‘ C25] INTZ:I<I+1 
[26] We x (Y-XB)x*(71xXI0))*x0.5 
C27] U0<((N,N)p1)xWxx(~1xXI0) 
[28] U1<+X1ixWxx("41xxI0) 
[29] U2+X2xyxx*(~1ixXxI0) 
C30] UU+R(3,MxN)p(,U0),(,U1),(,02) 
[34] JdJ0<U1<U2<10 
[32] Z<&, (WXXBxx(71xXI0) )+( ("1 xDELLA )+W x2 4H 
[33] W<10 
{34] B<ZhUU 
[35] 2<«UU<i0 
[36] BK+<B 
[37] XB8<«BC1]+(BL2]*X1)+B(3]xx2 
[38] RR<(Y-XB)x*x(71xxXI0) 
C39] XI0<R(N,M)pXI<XI+DIFFXI<(C-+/RRx ( (~1xDELTLA )+*RR) )+ 
((71xC)-4+/ (RR*2 )x*RR ) 
C40] +(L<2)/INTZ 
P C41] anaananaanaAnn VARIABILITY CHECK FOR UNCENSORED TIMES agnanaananan 
[42] VAR](+/+/ (71+ (Ut (&XB) )x (k71xXI0) x2 )t (T 14x) 
[43] U<i0 
a C4u] EPTA0<«~99999 
C45} X<Mpi 




























































J<0 . 
*(VARS1)/END 


- AMARAARAAAAAA INITIAL CONDITION FOR GAMMA PARAMETER ,ETA aannaana 


RRe (Y+XB.)%&* (7 1xXI0) 
 M1«70.5772 

M2<(+/+/DELTAXRR*2 )+ (4+/C) 

ETA0<0 

+(0>7+(( (ot )%*2)476 )+(M2%2 )=(M1x*2))/READY 
ETA0<" 1xe7 
READY :BADXI<0 

AA LARGE RECURSIVE LOOP INCLUDING ESTIMATE OF GAMMA PARAMETER ,ETA 
RECURSIVE: K+0 

J<dt1 

DATA+(S<+/(2]*RR),C 


FM*«ETAO FVALUE DATA 

AnAARAR BxSECTION SEARCH METHOD APPLIED. TO ESTIMATE ETA pananana 
BOUND: K+K+1 

FL«(LB+ETAO-Kx0.5) FVALUE DATA 

FR (RB*ETAO+KxX0.5)- FVALUE DATA 

AAAARAAAAAAARAAAAAR DETECT ILL CONDITION FOR ETA arnaaraAaAAAAAAAAAA 
+(K<10)/CHECKIRB 

+(v((( | FRO-FR)#0.5)<0.01)A(([FR)<0.1))/CHECKLRB 

ETA0+ 100000 

+END 
ARRARARAARAAARARAARARRARAAAAARAAARAAARRARARARAARARRARARARAAARARAR 


| CHECKERB: 


>(((xFM)x (XFL) )<0)/SELRB 

>(((*FM)x (xFR))<0)/SETLB 

FLO<FL 

FRO<FR 

>BOUND 
SETRB: RB<ETAO 

>BISECL 
SETLB: LB<ETAO 
BISECT: 

Fil<(MP+(LB+RB)+2) FVALUE DATA 

FL<+LB FVALUE DATA 
>((RB-LB)<(|0.01xMP))/REGRESS 
>((C(xFM)x(xFL) <0), CC (XFM)x (XFL) )>0))/SWAPR ,SWAPL 
SWNAPR: RB<UP 

>BISECL 
SWNAPL:3LB<MP 

*+BISECTL 

AnARA REGRESSION PROCEDURE FOR WEIBULL SCALE PARAMETERS ,BETA'S aa 
REGRESS: 

ETAO<(RB+LB)+2 
W<((CS<Q(N,M)p C(C+XETAO )¢S+xETAO ))xxRR)*0.5 
U0+((M,N)p1)xWx*(~1xXI0) 

U1<X1xWxx (~1xxXI0) 

U2*X2xWxx (~1xXI0) 

UU<& (3 ,MxXN)p(,U0),(,01),¢(,U2) 
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[97] 
~ £98]. 


[1007 


C102] 
[103] 


[99] 


[101] 


[104] 
[105] 


[106] 


C107] 
[108] 
[109] 


£110] 
[111] 


C112] 


[4113]: 


C144) 
[115] 
C116] 


(117), 


[118] 
[119] 
[120] 





Uoeu1+U210 


Z¢¥ , (WXXBx (~1xXI0) + (("1XDELTA ) 4x2 4H 
BeZQUU 

UU<Z<10: 
amanaaAAA NEWTON PROCEDURE FOR WEIBULL SHAPE PARAMETER,XI'S aan 
XB<BOL1]+(BL21xX1)+B031xxX2 

RR€(Y-XB)xx(71xXI0) 

Re (74x )+(+/02] RRx(("1xDELTA )+CSX*RR)) 

DR«(~1xC)=(+/02] CSx(xRR)xRRx2) 
nmaarAnAAAAAAAAAA DETECT THE ILL CONDITION OF XI'S anaananaaaana 
XX<e~(({H)<0.1)A((|DH)<0.01) 

>(~((x/XX )=0 )JABADXI=0 )/NEWLONXI 

BADXI#1 

X*#XX 
NEWLONXT 

XI0€8 (NM) pXI<XI+DIFFXI€¢XXx (~1xH DH 
an LEST FOR STOPPING CRITERIA AND UPDATE THE VALUE OF VARIABLES 
>((x/ (| (DIFPXI+#XI), (BK-B)+#B)<0.01)v (7250) )/END 

BK<B 

RR (Y-XB)xx(~1xXI0) 

*+RECURSIVE 
AAAARARAAA OUTPUT UNCENSORED LEVEL AND ALL ESTIMATES aagaaanaaan 
END: 

Re ((+/C)4MxN) , ETA0,8C1],802],B03], (C+/XIxX + (+/X)), +/0X 
V 
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